Dynamical mean-field theory of a double-exchange model with diagonal disorder 
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We present a simplified model for the colossal magnetoresistance in doped manganites by exactly 
solving a double-exchange model (with Ising-like local spins) and quenched binary disorder within 
dynamical mean field theory. We examine the magnetic properties and the electrical and thermal 
transport. Our solution illustrates three different physical regimes: (i) a weak-disorder regime, where 
the system acts like a renormalized double-exchange system (which is insufficient to describe the 
behavior in the manganites); (ii) a strong-disorder regime, where the system is described by strong- 
coupling physics about an insulating phase (which is the most favorable for large magnetoresistance) ; 
and (iii) a transition region of moderate disorder, where both double-exchange and strong-coupling 
effects are important. We use the thermopower as a stringent test for the applicability of this model 
to the manganites and find that the model is unable to properly account for the sign change of the 
thermopower seen in experiment. 
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I. INTRODUCTION 

An enormous interest in the family of doped man- 
ganese oxides Lai_a;Aa;Mn03 (in which A stands for Ca, 
Sr or Pb) has been created by the colossal magnetore- 
sistance (CM|R) exhibited in samples with doping levels 
in the rangeEJ 0.15 < x < 0.4. In such a doping region, 
there is a characteristic temperature Tp where the resis- 
tivity p has a peak; the CMR materials display metallic 
behavior (defined by dp/dT > 0) for T <Tp, and insulat- 
ing behavior (defined by dp/dT < 0) for T > Tp (except 
for LSMO at a; w 0.3). Hence there is a metal-insulator 
transition (MIT) or crossover at Tp. When placed in 
an external magnetic field H, the resistivity is strongly 
suppressed and the temperature at which the resistivity 
has a peak (Tp) increases. The magnitude of the relative 
magnetoresistance can be extremely large and can attain 
99% or more in some samples. 

It is widely accepted that the transport properties of 
these systems are closely related to their magnetic prop- 
erties. The temperature of the MIT is close to the Curie 
temperature Tc so one can say that these materials are 
metallic in the ferromagnetic phase and are insulating in 
the paramagnetic phase. The itinerant-electron and the 
local-spin ^^es are correlated by the double-exchange 
mechanismcTEl (which is one type of indirect exchange 
interaction between local spins via itinerant electrons). 
Double exchange consists of a cooperative effect where 
the motion of an itinerant electron favors the formation 
of ferromagnetic order of the local spins and, vice versa, 
the presence of ferromagnetic order facilitates the mo- 
tion of the itinerant electrons. Hence, the resistivity of a 
double-exchange system will increase when the temper- 
ature is increased through the Curie point. This is in 
qualitative agreement with the experimental data on the 
manganites. 

Unfortunately, it is also well known that double ex- 
change alone cannot explain the quantitative features of 



the temperature dependence of the resistivity through 
thej-entire temperature range observed in the mangan- 
itea3 (especially near the MIT). There are two proposed 
theoretical resolutions of this problem. The first one is 
that a large Jahn- Teller lattice distortionJs responsible 
for the anomalous transport propertieaJTij. The lattice 
distortion causes a metal-insulator transition via a strong 
(exponential) polaronic narrowing of the conduction- 
electron band. This polaronic narrowing also leads to 
a decrease of the Curie temperature because the double- 
exchange mechanism for ferromagnetic order is reduced 
when the bandwidth of the conduction electrons is nar- 
rowed. 

The polaronic mechaiusm has been criticized by a num- 
ber of authors. VarmatJ noticed that there exist double- 
exchange systems, such as TmSe^Tei-x, in which the 
transport anomalies seen in the manganites also occur, 
but a Jalm- Teller distortion is forbidden by symmetry. 
Furukawallj showed that a small-polaron picture leads to 
a strong suppression of the ferromagnetic transition tem- 
perature and the estimate of Tc due to double exchange 
was incorrectly calculated in Ref. ^ (Ref. |l^ reaches a sim- 
ilar conclusion). Further difficulties arise from the fajet 
that LSMO does not have a metal-insulator transitiontil 
at a; '--^ 0.3. Furukawa claims that LSMO [which has a 
relatively high value for Tc (about 380°K)] is a canonical 
double-exchange system (note that LSMO does showO a 
doping crossover from metallic behavior at x ^ 0.3 to ivi- 
sulating behavior at x ^ 0.15). Finally, an analysialj'EEl 
of the longitudinal and Hall resistivities in LCMO and 
LPMO cannot be explained in the small-polaron picture 
for the temperature range near Tp or for high tempera- 
tures (T»Tp). 

In the second proposed resolutiorO'Ej Ej, the insulat- 
ing behavior is caused by a combination of both magnetic 
disorder (due to the lack of ferromagnetic alignment in 
the paramagnetic phase) and nonmagnetic ionic disorder 
(due to the doping of the "A" metal) . The magnetic dis- 
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order arises from the "random" double-exchange factor 
cos 6/2 in the electronic hopping (where 9 is the angle 
between local spins). This is an offijdiagonal disorder that 
can lead to ^Xifchitz localizations of the charge carriers. 
Sheng et al.t§ll3 have found that this off-diagonal disor- 
der is insufficient to localize the electronic states at the 
Fermi level for moderate doping (which agrees with the 
claim that the double-exchange mechanism alone cannot 
describe the behavior of the manganites). The nonmag- 
netic disorder comes from the ionic doping of the 
ions (i.e. from randomness at the chemical substitution 
of La by A^+) which leads to a "random" local poten- 
tial for the charge carriers. This substitutional disor- 
der is always present in the doped materials, and it is 
physically meaningful to speak only about substitutional- 
disorder-averaged quantitiescJ. This ionic disorder is a 
diagonal disorder that can lead to an Anderson local- 
ization of tjhfi charge carriers. One-parameter scaling 
calculationalfl show that in the presence of a suitable 
strength of the ionic disorder, the magnetic disorder will 
cause the localization of electrons at the Fermi surface 
and induce a metal- insulator transi|tipn near T^.. (How- 
ever, there is experimental evidencec3 that Anderson lo- 
calization is not the cause of the metal-insulator transi- 
tion in Lao.erCas.saMnOs). 

Zhong et al.El have used dynamical mean field the- 
ory to study the metal-insulator transition in the man- 
ganites in the framework of an s — d model with clas- 
sical local spins and doping-induced disorder (see also 
Ref. ^3|). They were able to show that a MIT is possible 
for a binary-alloy distribution of the ionic energy lev- 
els due to a splitting of the electronic band (correlated 
by the double-exchange process) into completely filled 
and empty subbands at some critical value j-of disorder 
strength. Such a distribution was also usecta to study 
the Hubbard model with substitutional disorder within 
dynamical mean field theory. 

In this contribution, we consider the magnetic and 
transport properties of a simple double-exchange system 
with diagonal disorder. The system is described by the 
following Hamiltonian 



where the c-operators are composite operators 



-(1 + asl)ai„, 



(1) 



(2) 



with si the z-component of the local spin (5 = 1/2) 



described by the diagonal Pauli matrix, (s 



z\2 



1, and 



aj^(aj^) the ordinary Fermi annihilation (creation) oper- 
ator for an itinerant electron with spin projection a at lat- 
tice site i. The first term in Eq. (|^) describes the doping- 
induced (diagonal) disorder (the energies are chosen 
from a disorder distribution) and the second term rep- 
resents a simplified version of the quantum-mechanical 
double-exchange mechanism for ferromagnetic ordering 



of local spins. This term can be obtained from an s — d 
model with local spins that are described by Ising spins in 
the limit of infinitely strong exchange interaction between 
itinerant and localized electrons. This simplification of 
describing the local spins by Ising variables conserves 
the main feature of double exchange: the second term 
in Eq. (|^) only allows the transfer of itinerant electrons 
with spin parallel to the local spin at every site of the lat- 
tice [see Eq. (^]. But, these simplified operators do not 
allow any spin-flip processes. Such processes can be im- 
portant at low temperatures where the thermodynamics 
of the system is governed by spin- wave excitations. How- 
ever, since dynamical mean-field theory cannot describe 
spin waves, including such effects is beyond the scope of 
this work. Spatial correlations between the local spins of 
the form (sf s|) [which is containet in Eq. (|l])] should be 
important near the Curie pointErEZi Tc, but they are also 
beyond the capabilities of dynamical mean field theory. 

This paper is organized as follows: The dynamical 
mean-field theory equations for the system are presented 
in Section II. The binary probability distribution for the 
doping-induced disorder and simplifications for this dis- 
order averaging are discussed in Section HI. In Section 
IV, the influence of disorder on the magnetic properties 
(decreases of Tc, the paramagnetic susceptibility, and 
the magnetization of the local-spin subsystem) are in- 
vestigated. Here, we show that for strong disorder, the 
double-exchange mechanism of ferromagnetic ordering is 
replaced by a disorder-induced ferromagnetism and we 
describe three characteristic disorder regimes: (i) the 
weak-disorder regime where the double-exchange mech- 
anism dominates; (ii) the strong-disorder regime where 
the magnetic and transport properties are determined 
by strong-coupling physics from an insulating phase; 
and (iii) a transition regime where both mechanisms are 
important. The main transport properties (resistivity, 
thcrmopower, and thermal conductivity) in these three 
regimes (along with the magnetoresistivity) are discussed 
in Sections V and VI, respectively. Section VII contains 
our concluding remarks. 



II. FORMALISM FOR THE DYNAMICAL MEAN 
FIELD THEORY 

The dynamical mean field theory equations for the sys- 
tem described by Eq. (|^) can be obtain io-two ways: (i) 
a diagrammatic technique for|-ei-operatorsEj can be com- 
bined with disorder averaging23 in the limit z ^ oo (z is 
the coordination number), or (ii) one can work directly 
with a local effective action Sef / . The first method yields 
a direct calculation of the dynamical mean field equations 
for the disorder-averaged band Green's function and for 
the magnetization of the local-spin subsystem by exactly 
summing all nonvanishing graphs as z — > oo. This pro- 
cedure is cumbersome, so we will use the effective action 
approach here. 
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Since the anticommutator of c-operators is not a num- 
ber, it is convenient to begin with the original s — d model 
with Ising spins and diagonal disorder: 



i i(7 



(3) 



where we have introduced two external magnetic fields 
{h acts on the local-spin subsystem and H acts on the 
itinerant-electron subsystem). This is done to allow 
derivatives with respect to the fields to be calculated 
properly. In the end results, we set the two fields equal to 
the true external magnetic field. We will take the limit 
where the s — d exchange parameter / becomes infinitely 
large (/ cxd), because this is the regime where the 
s ~ d Hamiltonian is mapped onto the double-exchange 
Hamiltonian. 

The mathematical structure of the s — d Hamiltonian m 
Eq. is similar to that of the Falicov-Kimball modeled 
which can be solved exactly in infinite dimensionsEil. The 
procedure is to first solve the atomic problem in an ex- 
ternal time-dependent field and then adjust the field so 
that the atomic Green's function equals the local lat- 
tice Green's function. The local effective action for this 
atomic problem is 

Seffie) = -\l3hs^ draUr) [e-f,- ^aH 



- iw^a,(r)-f ^ ^dr ^ dr' a\(T) 

'^0 

X (^|:5(r-r') + A.(r-r'))a,(r'), (4) 

with /3 = l/T and Ao.(t) is the time-dependent field. The 
partition function becomes 



(5) 



where the trace is taken over s^. The disorder-averaged 
free energy becomes 



Fe// = ^T(lnZ,;/(e)) ^ -T ^ rfeF(e) lnZe//(6). 



(6) 



Here, P(e) is a probability distribution function for the 
(random) atomic energies and the angle brackets (. . .) 
denote the disorder averaging. 

The disorder-averaged local band Green's function 
(Cfna) is determined by a functional derivative with re- 
spect to the atomic field. In Fourier space, we have 



dF, 



eff 



dAno 



(7) 



where A„a = Aa-{iujn) is the Fourier transform of A^{t) 
and cun — (2n -|- l)7rT is the Fermionic Matsubara fre- 
quency. 

In taking the limit / — )■ cx) in Eq. (^, we must first 
renormalizc the chemical potential (/i -I- //2 — > ^). Then 
Eq. (0) becomes 



{Gn 



deP{e)' 



^[1 



(8) 



where 



a„o- = it^n + M - A„cr, (9) 

is the inverse of the effective medium and 

m{e) =t'Anh^[f3h + f3H/2 + Xpie)], (10) 

is the magnetization of the local spins (when the band 
electrons have disorder energy e) and 



The total magnetization is defined to be 



eff 



d{h/2) 



= j deP{e)m{e). 



(11) 



(12) 



Note that the expression in Eq. (|T^) for the magnetiza- 
tion contains a hyperbolic tangent, just like in the mean- 
field theory of an Ising ferromagnet. Therefore, Ai?(e) 
can be identified as an internal molecular field acting on 
the local spins. Eqs. ( pi] ) and (||) show that this molec- 
ular field is completely determined by the properties of 
the itinerant-electron subsystem. Hence, the magnetic 
and transport properties of the double-exchange system 
are correlated. 

In infinite dimensions, the inverse of the effective 
medium ana, local self energy E„ct = So-(iWn), and the 
local Green's function {Gna) are related by 



{Gn 



(13) 



The self consistency relation equates the atomic Green's 
function with the local Green's function of the lattice. 
The latter can be calculated from the local self energy 
by summing over all momentum. Since the self energy is 
independent of momentum, we find 



{Gna) ^ JjYl " / 

7, 



dxD°{x) 



ioJn -I- /X - S„ 



(14) 



where D° {x) is the density of states of the noninteracting 
itinerant electrons. We choose to examine the infinite- 
coordination Bethe lattice, where 
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(15) 



with 2t* — 1 chosen to be the energy unit. The integral 
in Eq. can be computed exactly yielding 



*2 



2t* 



2t* 



(16) 



Replacing I]„cr by a„cr — (^Jncr) ^ and solving for the 
Green's function gives 



t*'^{Gna) = iw„ + /i - a„ 



and 



*2 



(17) 



(18) 



Substituting (17) into (g), we obtain the following final 
equation for ana- 



ILOn + fl ~ t* 



^[1 



deP{e)^^-^-f^. (19) 



This equation is simpler than the corresponding system 
of equations analyzed by Zhong et alE2l. However, we 
believe that our approach captures the main features of 
this system, as we described in Section I. 

The chemical potential ^ is adjusted to give the correct 
filling for the itinerant electrons: 

1 f°° 

n = ^ / dujf{uj)lra{ga {uj - ji)) 

oo 

= ^J2j dLjf{Lo)lmaa{u;-fi), (20) 



where 



1 



exp /3(ijj — ii) + 1 



(21) 



is the Fermi-Dirac function and aa-{uj) is the solution of 
Eq. ( p^ ) (evaluated on the real- frequency axis). 

Using Eq. (p^), we can also evaluate both the magne- 
tization m and the paramagnetic susceptibility x of the 
local-spin subsystem: 



dm 



h=H=m=0 



deP{e) 



dm{e) 



dh 



(22) 



h=H=m=0 



where dm{e)/dh is defined by the following integral equa- 
tion 



dm(e) 
dh 



9^ 9 ^ n _ 



1 



, F(e') dm(e') 



with 



1 ~ An J dn — e' dh 



An - It*' I &P(e)^^-^ 
2 J (an - ey 



(23) 



(24) 



We note that an integral equation for a two-particle cor- 
relation function is expected for disO|Cdered systems de- 
scribed by the Falicov-Kimball modelo. 



III. BINARY PROBABILITY DISTRIBUTION 

An analysis-ofn the electronic properties of the man- 
ganites show£^'E3 that a binary distribution for doping- 
induced disorder can be approximately suitable for the 
doped materiaJs, This distribution is written in a sym- 
metric form asL^I 

P{ei) - (1 - x)5{e, + ^A) + x5{e, - ^A), (25) 

where x is the fraction of the sites having an additional 
local potential A (disorder strength) due to the ionic 
doping [the symmetric form requires a renormalization 
of the chemical potential ^ ^ n — {1 — x)/S./2\. It is 
seen from Eq. (|2^) that this choice for the distribution 
behaves electronicallyj-as-a coherent superposition of its 
end-point compoundsEllH. 

If the chemical substitution of the ions also causes the 
appearance of holes in the band, then there is a correla- 
tion between the electron density and the concentration 
of dopant sites 



1 



n 



(26) 



where v is the number of electronic states per lattice 
site. In the system we consider, double occupation by 
itinerant electrons is excluded so ^' = 1. Therefore, we 
have 1 — a; 71 in the probability distribution of Eq. (|2F 
This constraint is important to allow a MIT. 

The basic dynamical-mean-field equations are simpli- 
fied for the binary distribution. In particular, Eq. 
becomes 

1 *2f(l -")[!+ ^™(iA)] 



n[l + o-m(-iA)] 



iai/-iA 



(27) 



where m(A/2) and m(— A/2) depend on the complete 
set {ana}- Aside from the factors of [1 -I- crm(±A/2)]/2, 
this result is identical to-that of an annealed binary alloy 
problem, as first solvedEil by Brandt and Mielsch. 

The main feature of the probability distribution ( [25| ) 
[with \ ^ X = n] is that the chemical potential is located 
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in the gap for any electron density (in the paramagnetic 
phase) when the conduction band is spUt by strong disor- 
der A > Ac. Indeed, when A = 0, the local band Green's 
function has the form 

(28) 

This shows that the bandwidt h of the p ure double- 
exchange system is equal to At* -^(1 -I- to)/2 for the (ma- 
jority) spin-up electrons and 4i*-\/ (1 — m)/2 for the (mi- 
nority) spin-down electrons. At T = (where m — 1), 
there are no spin-down electrons (the spin-down electron 
bandwidth is zero), and the spin- up electrons act like free 
electrons with the full value of At* — W for the band- 
width. This ferromagnetic ordering promotes the motion 
of the electrons. Increasing the temperature destroys the 
ferromagnetism; the bandwidth for the spin-up electrons 
is decreased and the bandwidth for spin-down electrons is 
increased so that they are both equal to At* \/l/2 in the 
high-temperature paramagnetic phaseE3. In the param- 
agnetic phase the bandwidth, and density of states are 
independent of temperatureEJ. 
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FIG. 1. Density of states as a function of frequency for 
different disorder strengths A (in the paramagnetic phase). 
Vertical lines indicate the location of the chemical potential. 
Note that the density of states is plotted on an absolute en- 
ergy scale, so the chemical potential is not shifted to lie at 
ui = 0. The inset displays the dependence of the critical val- 
ues of disorder in the paramag netic [A(r = T^) = Af ] and 
ferromagnetic [Ac|(T = 0) = Af ] phases as a function of the 
electron filling n. 

Figure 1 shows the influence of disorder on the 
conduction-electron density of states in the paramagnetic 
phase at n = 0.8. When the strength of the disorder is 
small, the band is only slightly distorted from the semi- 
circular shape. As the disorder increases, a pseudogap 
first appears and then a true gap develops when the dis- 
order is larger than a critical value. Because the double- 
exchange couples ferromagnetism to the mobility of the 
itinerant electrons, the critical value of the disorder de- 
pends on the spin polarization of the system and on the 



total electron density. This is depicted in the inset to Fig- 
ure 1, where the critical value of the disorder, required 
for the metal-insulator transition, is plotted as a function 
of the total electron density. Two curves are shown: (i) 
the critical disorder strength (Af ) needed for the transi- 
tion in the high-temperature paramagnetic phase T > T^, 
where there is no spin polarization (bottom curve); and 
(ii) the critical disorder strength (Af ) in the fully po- 
larized ferromagnetic phase at T = (top curve). Since 
the ferromagnetic order makes the bandwidth larger, the 
critical value of disorder needed for the MIT increases 
as one enters the ferromagnetic phase. If the disorder 
is large enough, it is always an insulator (even in the 
ferromagnetic phase) — but there is a regime, where the 
system can be an insulator in the paramagnetic phase, 
and a metal in the ferromagnetic phase (at least it is 
metallic for the majority [spin-up] electrons, it would be 
insulating for the minority [spin-down] electrons). This 
is the regime that is relevant for the CMR materials. 

In the limit where A ^ oo, the band is always split, 
and the upper subband is pushed to infinite energy and 
can be neglected. From Eq. (^^ we have 

Ima(w - ii) = iv/2t*2n - w^, < 2t*'^n, (29) 

for H = {). The magnetization vanishes, because the 
internal molecular field \p vanishes when A — > oo as 
seen from Eq. (pT[). Hence, there is no ferromagnetic 
order, and the inverse of the effective medium satisfies 
a^iuj) — ai{uj) ~ a(uj). Since 

-J^ j dujln\a{uj - /i) = n, (30) 

for this case, the lower band is completely occupied and 
the chemical potential lies in the gap for all n. 

For finite A, the same conclusion is obtained for the 
high-temperature paramagnetic phase. This follows di- 
rectly from the numerical calculations, but can be under- 
stood from the fact that the total spectral weight in the 
lower band doesn't change until the gap closes. The fer- 
romagnetic transition temperature generically increases 
from zero though. Figure 1 shows the location of the 
chemical potential at Tc with vertical lines for A = 0, 
A = 0.4 and A = 0.6 (note that in the last case the 
chemical potential lies in the pseudo-gap). For A = 0.8 
and A = 1.0 (Ac = 0.6581 at n = 0.8), the chemical 
potential lies in the gap. 

This situation is similar to what takes place in the 
static Hols tcip .model, which is solved exactly in infinite 
dimensionsE3£3. In this model, bipolaron formation leads 
to the opening of a gap when the electron-phonon inter- 
action is strong enough, and the chemical potential lies 
in the gap for all n. It differs from the Mott-Hubbard 
transition which occurs only at half-filling. 
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IV. MAGNETIC PROPERTIES 

The magnetism of the binary-disorder double-exchange 
system is determined by the uniform (ferromagnetic) 
magnetization 



m = {1 - n)m(^^A^ + nm(-^Ay (31) 
and by the uniform (ferromagnetic) susceptibiUty 



.cM|A) dm(-iA) 
X = (1 — n) rr h 



dh 



dh 



(32) 



The algebraic equations for dm{A/2)/dh 
and dm{~A/2)/dh are taken atft. = iJ = m = Oin 
Eq. ( p3| ) with the binary probability distribution for the 
disorder. The determinant of this system of equations 
equals zero at the ferromagnetic Curie temperature Te- 
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FIG. 2. Curie temperature as a function of disorder 
strength A and electron filling n: (a) shows results for con- 
stant electron filling, while (b) shows results for constant dis- 
order strength. 



Figure 2 shows (a) the disorder dependence of the 
Curie temperature for different electron densities and (b) 
the electron density dependence of the Curie temperature 
for different disorder strengths. The value of Tc for the 
pure double-exchange system [the curve A = in Fig- 
ure 2(b)] essentially coincides with that obtained by Eu,- 
rukawa for double exchange with classical local spinsll3. 
One needs to choose the electronic bandwidth to be quite 
narrow {W on the order of 0.5 — 1 ev for Tc ranging 
from 315 — 630°K at n ~ 0.5). These estimates of Tc for 
the pure double-exchange system are comparable with 
some materials (like LSMO where Tc ~ 380°K), and are 
much smaller than those predicted in Ref. |^. Neverthe- 
less, if we take a more reasonable value for the bandwidth 
(W ~ 2 ev), then it is clear that double exchange alone 
cannot explain the values of Tc for the manganites. 

Figure 2 shows that disorder suppresses the ferromag- 
netic transition temperature. The physics of this is clear. 
Carrier motion promotes the double-exchange ferromag- 
netic order, whereas disorder reduces the electron mo- 
tion, and thereby it reduces the ability of the double- 
exchange process to produce ferromagnetism — the net ef- 
fect is to reduce Tc- The calculated disorder dependence 
of Tc agrees qualitatively with that found for the com- 
bined double-exchange-Holsteip-jnodelJ and with Nari- 
manov and Varma's calculatiorcS with a Gaussian disor- 
der probability distribution. 

Figure 2 shows also that three disorder regimes can be 
distinguished for our model. The first regime (renormal- 
ized double-exchange) is the regime where Tc depends 
only weakly on disorder and corresponds to the flat re- 
gions of the curves in Figure 2(a) [A < 0.5]. We find that 
Tc is suppressed most strongly at 0.5 < A < 1.0 [see Fig- 
ure 2(a)]; i.e. in the vicinity of the critical values of disor- 
der where a gap in the density of states is created. We call 
this regime the transition regime, where the ferromag- 
netism process is changing from a renormalized double- 
exchange process to a renormalized strong-coupling pro- 
cess. In this range of A, two crossovers occur: (i) the 
metallic conductivity is replaced by a thermally activated 
(insulating) conductivity, and (ii) double-exchange ferro- 
magnetism is replaced by a ferromagnetism that is caused 
by virtual electron transfers from the filled lower band to 
the empty upper subband (and vice versa). Indeed, the 
double-exchange mechanism for the ferromagnetic order- 
ing of the local spins is caused by real electron transfers, 
which are possible only at small disorder. For strong 
disorder (A > 1), when the itinerant electrons are lo- 
calized, the ferromagnetic ordering is caused by virtual 
transfers to neighboring sites (that have an additional 
local potential A) and back. (All sites without this ad- 
ditional potential are filled so hopping onto those sites is 
impossible.) Hence, the Curie temperature is inversely 
proportional to A at strong disorder [the 1/A behavior 
is clearly seen in Figure 2(a) for A > 1.5]. 

Usually, virtual electron transfers lead to a Heisenberg 
type of ferromagnetic ordering, i.e. the magnetic en- 
ergy contains a cos & dependence, not the cos Q /2 depen- 
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dence of double exchange. Therefore, the crossover from 
double-exchange ferromagnetism to disorder-induced fer- 
romagnetism (with a Heisenberg type of magnetic en- 
ergy) must be seen in the changing of the behavior of 
magnetic quantities (for example, paramagnetic suscep- 
tibility and magnetization) as the disorder strength in- 
creases. 

In particular, an increase of the disorder strength 
changes the temperature dependence of the uniform sus- 
ceptibility. Indeed, the pure double-exchange system 
(A = 0) reveals Curie-type behavior for x at high tem- 
peratures (x ~ l/T')) i-e. the Curie temperature is equal 
to zero. As the temperature is decreased, the curvature 
of the inverse susceptibility changes, so that it inter- 
sects the temperature axis at a finite temperature, yield- 
ing a nonzero value for and a Curie- Weiss law for x- 
This feature of double exchange was noted in the pioneer- 
ing work by Anderson and HasegawaB (see also Ref. 28). 
(It should be noted that this aspect of double exchange is 
still not completely understood. See, for example, Ref. ^ 
where a high-temperature expansion is employed). 

Numerical calculation of the uniform susceptibility, 
Eq. (^^, at strong disorder shows that x obeys a Curie- 
Weiss law with a nonzero Curie temperature in this case. 
Thus, the crossover to a disordered-induced mechanism 
of ferromagnetic ordering leads to a change in the tem- 
perature dependence of x- 
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0.020 



FIG. 3. Magnetization of the local-spins m as a func- 
tion of the external magnetic field H for difi'erent disorder 
strengths A. The inset shows the disorder dependence of m 
aX H = 0.005 and H = 0.01 for the double-exchange model 
(solid lines) and for the Ising model in the mean-field approx- 
imation (dotted lines). 

In Figure 3 we plot the magnetization m as a function 
of external magnetic field H at T = Tc for a number 
of different disorder strengths. The field-induced mag- 



netization at fixed magnetic field depends on disorder 
and the inset to Figure 3 shows that the magnetization 
(solid lines) initially decreases as the disorder strength 
increases until one reaches the region where the conduc- 
tion band has a well-developed pseudogap and begins to 
split (Ac = 0.6581 for n — 0.8). As the disorder increases 
further, the magnetization starts to increase and it can 
be described by the following mean-field equation 



m = tanh/5( mT^ + 



(33) 



for the magnetization of an Ising model (or a Heisenberg 
model with S = 1/2). Dotted lines in the inset show the 
field-induced magnetization evaluated from ( |33| ) where 
Tc is equal to the Curie temperature of the disordered 
double-exchange model at a given A. Hence, the disorder 
dependence of the magnetization at fixed magnetic field 
also shows the crossover from double-exchange-induced 
to disorder-induced ferromagnetism at strong disorder. 




0.000 0.005 0.010 0.015 0.020 

IVIagnetic field, H 

FIG. 4. Magnetization of the local-spins, m, as a function 
of the external magnetic field H for different relative temper- 
atures T = T/Tc. 1 - r = 0.4, 2 - T = 0.8, 3 - r = 1.0, 4 - 
T = 1.1, 5 - r = 1.2, 6 - r = 1.3, 7 - r = 1.4. 

From Figure 3 we see that the most favorable condi- 
tions for a large field-induced magnetization are strong 
disorder and strong magnetic fields. Figure 4 shows the 
field dependence of ni for different relative temperatures 
r = T /Tc and fixed disorder. The largest absolute value 
for m is observed at small t {t — 0.4) but the relative 
increase of m (when H is also increased) is weak. The 
strongest growth of m at weak fields is seen near the 
Curie temperature. At T — Tc {t = 1.0) this growth is 
maximal, but even at t — 1.4 the H dependence of m 
is linear. Note that such behavior of the field- induced 
magnetization is in good aepeement with the experimen- 
tal data on the manganitegij. 
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Our analysis of the magnetic properties has shown 
the emergence of three different regimes: (i) a renormal- 
ized double-exchange regime (or weak-disorder regime); 
(ii) a strong-coupling physics regime (or strong-disorder 
regime); and (iii) a transition regime where the proper- 
ties crossover from metal to insulator and from double- 
exchange ferromagnetic ordering to strong-coupling- 
induced ferromagnetism. This classification scheme will 
be more sharply defined as we examine the transport 
properties in Section V. 



V. TRANSPORT PROPERTIES 

In this section, we consider the electrical and thermal 
transport properties of the double-exchange system in 
zero external magnetic field. The transport is determined 
from the following set of equations^. The electrical con- 
ductivity cr = 1/p is 



a = Lii 



the thermopower is 

ks 1 Li2 



\e\ TLn' 



(34) 



(35) 



and the thermal conductivity (of the electronic system) 
satisfies 



e 



where the transport coefficients Lij are defined by 

oo 

J d£t;2 (£)£)"(£) 



Li. = na-Q 



duj 



{u;-tiy+^'^Al{e,u;-ti). (37) 



Here the spectral function is given by 

1. 



(38) 



(To is the unit of conductivity, t/(j(e, w) is given in Eq. (|l4 
and v{e) is the current vertex which is equal to 



V(4i*2-£2)/3^ (39) 

for the Bethe latticeSS with z ^ oo. 

Figure 5 shows the temperature dependence of the re- 
sistivity for different electron fillings in the three differ- 
ent disorder regimes. (In this and the following Figures, 
the resistivity is plotted in units of po = l/uo)- In Fig- 
ure 5(a), the weak-disorder regime, the transport proper- 
ties are determined by the double-exchange mechanism. 



From Eqs. (y_8|) and (^Tj), one obtains the following ex- 
pression for the electronic self-energy 

1 / 1 — am 



2 \ 1 + (TTO^ 

X |t^-HAi± V(^^ + M)^-2t*2(l + am)}, (40) 
at A = 0. Therefore, 



iial^a- [uj + id) = 



(41) 



for {uj + /i)2 < 2^*2(1 + am) and lmJ:„{uj + iS) = Q 
otherwise. 

Note that the factor 



1 — am 
1 -|- am ' 



(42) 



in Eq. ( pl| ) is contained in the expression for the self- 
energy in Ref. |4^ for the pure double-exchange system 
with classical local spins and a Lorentzian density of 
states (we use a semicircular density of states here). 
This factor plays an essential role in the low-temperature 
behavior of the resistivity [Figure 5(a)]. It provides a 
decrease in the resistivity when the magnetization in- 
creases. Indeed, in the limit m ^ 1 or T — > 0, we have 



(43) 



and the conduction-electron subsystem becomes a free- 
electron gas of spin-up electrons with a conductivity that 
is a delta function at zero frequency (whose strength is 
the Drude weight). 

On the other hand, in the high-temperature param- 
agnetic phase, m = 0, and the correlated bandwidth is 
narrowed by a factor of \/2 due to paramagnetic spin 
disorder. In this case, the expression for the self-energy 
exactly coincides with the one obtained in Ref. ^ for the 
pure double-exchange system. Since the density of states 
is independent of temperature here, all of the tempera- 
ture dependence of the resistivity arises from the Fermi 
factor [—df{uj)/diL;] [see Eq. (p7|)]. If one makes the as- 
sumption that the chemical potential is also temperature- 
independent, so that the derivative of the Fermi factor 
depends . ouW weakly on temperature, then one would 
concludeOtS that the resistivity in the paramagnetic 
phase (due to double exchange only) is essentially a con- 
stant. However, more accurate calculations, that take 
into account the temperature dependence of the chem- 
ical potential, show that the pure double-exchange sys- 
tem becomes a bad-metal with dp/dT > [see, for ex- 
ample, Ref. 1^ where the bad-metal behavior was shown 
in a double-exchange system with weak electron-phonon 
interaction in their Figure 6(a)]. 

The sharp decrease of the resistivity at T < Tc for all 
n is caused by the rapid increase of the magnetization 
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m. This discontinuity in the slope dp/dT at T = Tc 
is a consequence of the dynamical mean-field approach. 
Incorporation of spatial spin fluctuationg23 will smooth 
out the temperature dependence of the resistivity in the 
vicinity of T^, but this is beyond dynamical mean field 
theory. 
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FIG. 5. Temperature dependence of the resistivity p for 
different disorder strengths A: (a) A = 0, (b) A = 0.5, (c) 
A = 0.75 and (d) A = 1.0. The curves labeled by 1 corre- 
spond to an electron filling of n = 0.6; 2 denotes n = 0.7; and 
3 denotes n = 0.8. 

As disorder is added to the system, the properties are 
initially changed little. When the disorder becomes large 
enough A « 0.4, then we start to feel a more direct in- 
fluence of the disorder as it produces a pseudogap in the 
density of states and we enter the transition regime of 
moderate disorder. In this regime, the slope of the resis- 
tivity can become negative dp/dT above Tc as shown in 
Figures 5(b) and 5(c). The value A ~ 0.4 is the boundary 
value that separates the weak-disorder and the moderate- 
disorder regimes for n = 0.8. As the disorder is in- 
creased further, the metallic conductivity is gradually re- 
placed by a thermally activated conductivity (this starts 
at A « A^). As one enters the strong-disorder regime, 
there is a MIT at Tc. Here the system displays insulating 
behavior everywhere, except just below the Curie point, 
where the rapid increase in the magnetization can cause 
the resistance to drop over a small temperature range 



before it turns around and increases again. This occurs 
in the transition from the paramagnetic insulator to the 
ferromagnetic insulator because the charge gap in the 
ferromagnetic insulator is smaller than the charge gap in 
the paramagnetic insulator. Those sharp cusps seen in 
Figure 5(d) will generically be smoothed out by spatial 
fluctuations. 




0.2 0.4 0.6 

Disorder strength, A 

FIG. 6. Resistivity at T = (dotted line) and at T = 
(solid line) versus disorder strength. Inset is a plot of 
the Curie temperature Tc versus the residual resistivity 
Pr = p{T = 0). Note how the two resistivity curves cross at 
a critical value of disorder, which determines the beginning 
of the transition region from the moderate to strong-disorder 
regimes. 

We can analyze the transition from moderate to strong 
disorder more quantitatively. We focus on two character- 
istic resistivities: the resistivity at Tc and the resistivity 
at r = (residual resistivity). For weak disorder the 
residual resistivity is much smaller than the resistivity 
at Tc- When we reach the moderate disorder regime, 
the residual resistivity starts to increase rapidly, eventu- 
ally overtaking the resistivity at Tc in the strong-disorder 
regime. We denote the boundary between the moder- 
ate and strong-disorder regimes by the value of disorder 
where p{T = 0) = p{T = Tc). This occurs at A ~ 0.892 
for n — 0.8 as shown in Figure 6. This transition occurs 
when the density of states has developed a strong pseu- 
dogap, but has not yet become an insulator (Af = 0.931 
for n — 0.8). The range of disorder between these two 
limits, (0.892 < A < Af = 0.931 for n = 0.8) is the 
transitional region between the moderate-disorder and 
the strong-disorder regimes. 

This analysis is in qualitative agreement with tJie cal- 
culations that use a Lorentaian density of stateso and 
especially with calculations^ that show a divergence of 
the T — resistivity at a critical disorder strength. Note 
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that Anderson localization is the cause of the MIT in 
Ref. 
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FIG. 7. Temperature dependence of (a) the thermopower 
and (b) the electronic thermal conductivity for different values 
of disorder. 

The inset to Figure 6 shows the Curie temperature 
as a function of the residual resistivity p{T = 0). This 
plot was obtained by combining the A dependences of 
p(r ~ 0) and Tc (see Figure 2). It is clearly seen that 
the suppression of due to disorder is accompanied by 
the increase of the residual resistivity. We find that our 
functional dependence of Tc on p{T — 0) is smoother 
than tliat,found in Ref. Ilj and agrees better with exper- 
imentBEl. 

Now we examine the thermal properties of our system 
including the thermopower S'(T) and the electronic ther- 
mal conductivity k(T) which are shown in Figure 7. The 



behavior of these two quantities are quite different from 
each other — while the thermal conductivity always van- 
ishes as r — !■ 0, the thermopower will either vanish or 
diverge depending on whether the system is metallic or 
insulating as T ^ 0. The thermal conductivity behaves 
as expected with a sharp increase at Tc due to the opening 
of conduction channels as the magnetization grows, and 
a linear decrease to zero at low temperatures. Note that 
the thermal conductivity vanishes at low temperatures 
even in the metal, because the heat current vanishes. 

The thermopower behaves quite a bit differently. It too 
shows a strong effect at Tc (in this case a sharp decrease), 
but the low temperature behavior is most interesting. 
For weak disorder (metallic phases), the thermopower 
vanishes as T — > 0, but once a gap opens in the density 
of states, the thermopower diverges as T ^ since the 
chemical potential lies in the gap. Peltier's coefficient, 
P = TS(T), is approximately given by a straight line for 
strong disorder (A = 1.0) and low temperature (where 
the magnetization has saturated m ~ 1), i.e. P — a + 
7T, a ~ 0.17/iV. Hence, the thermopower S{T) can be 
approximately represented by 



(44) 



for low temperature and strong disorder. This relatioa 
is typical of what is seen in intrinsic semiconductorscSI 
Note that the thermopower behaves similarly forj-elec- 
tronic systems that undergo Anderson localizations: if 
the chemical potential is in the region of localized states, 
then S{T) also diverges as T ^ 0. 

The slope of the thermopower dS/dT has a disconti- 
nuity at Tc, hxit-S{T) does not change sign in our model. 
The prediction^ that S{T) alter its sign at T = 
in double-exchange systems was founded on two prin- 
ciples: (i) the itinerant-electron subsystem is a Fermi- 
liquid and (ii) the derivative of the chemical potential 
dp/dT changes sign at T = Tc. While we find that the 
derivative dfi/dT ^does change sign at T = Tc in agree- 
ment with otherstaES, the itinerant-electron subsystem 
is not a Fermi-liquid in our model, because the imagi- 
nary part of the self-energy does not vanish as T — ^ at 
the Fermi-energy. Hence, the reasoning that led to the 
prediction of a sign change cannot be applied here. 

The high-temperature behavior of the thermopower for 
strong disorde^is similar to that of S{T) in a small- 
polaron modelcSI where the conductivity is thermally acti- 
vated and S{T) ~ ln[c/(l — c)]; with c the small-polaron 
concentration and c < 1. Indeed, the thermopower has 
a weak temperature dependence for strong disorder [see 
Figure 8(a)] in the paramagnetic phase and T ^ Tc. Fur- 
thermore, for fixed temperature, we find that S{T) de- 
creases when the electron filling is decreased to n — 0.5. 
At n = 0.5 the thermopower is equal to zero at all tem- 
peratures as in the small-polaron theory. Below n = 0.5, 
we find that the thermopower changes sign. Hence we 
only find a sign change of the thermopower when the 
electron filling is varied. 
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VI. MAGNETORESISTANCE 



ture qualitai;iiifily agrees with the experimental data on 
manganiteada. 



The most interesting property of the manganites is 
the fact that the resistance changes so dramatically in 
a magnetic field. This makes them useful as possible 
magnetic field sensors for the magnetic storage commu- 
nity. We find similar magnetoresistance effects in our 
model, especially when we are close to the Curie point. 
The origin of the magnetoresistance lies in the sensitivity 
of the resistivity to the magnetization, and the ease with 
which m can be tuned by a magnetic field. Typically, 
the field increases m, which then reduces the resistivity. 
Experimentallyfij, there is a strong correlation between 
the field-induced changes in p and m. 

In order to calculate the magnetoresistance, we must 
first solve for the conductivity in the presence of an exter- 
nal magnetic field. If we perform a simple shift in the def- 
inition of the inverse effective medium a„o- — > ana+Ha/2 
and replace the integration variable in Eq. ( p7[ ) by w — > 
u; + Ha/2, then the only modification is that the deriva- 
tive of the Fermi factor [—df{uj)/dLu] now has an explicit 
H dependence and the spectral function A^- in Eq. ( |29| ) 
depends on H through the magnetization m. The direct 
dependence of p on if via the factor [—df{uj)/duj] always 
yields a small positive magnetoresistance that is caused 
by the Zeeman interaction with the external magnetic 
field; it can be neglected in weak fields. 




0.000 0.005 0.010 0.015 0.020 



Magnetic field, H 

FIG. 9. Resistivity p{T) as a function of external mag- 
netic field H for different relative temperatures r = T/Tc'. 
1- r = 0.4, 2 - r = 0.8, 3 - r = 1.0, 4 - r = 1.1, 5 - r = 1.2, 6 
- r = 1.3 and 7 - r = 1.4. 
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FIG. 8. Temperature dependence of the resistivity p{T) for 
different magnetic fields H: I - H = 0.0, 2 - H = 0.005, 3 - 
H = 0.01, A - H = 0.015, 5 - if = 0.02. 

Figure 8 shows the temperature dependence of the re- 
sistivity for different magnetic fields and two disorder 
regimes: (a) moderate-disorder and (b) strong-disorder. 
The field-induced modifications to the magnetization 
suppress the resistivity. This effect is strongest in the 
vicinity of the Curie point, because the spin susceptibil- 
ity is large there; hence the field changes the magneti- 
zation most strongly there. The peak position in p{T) 
shifts to higher temperature with increasing H and there 
is a critical value of H above of which the MIT disap- 
pears (i.e., dp/dT > for all temperatures). This pic- 



The magnetic-field dependence of the resistivity at dif- 
ferent relative temperatures t = T/Tc (Figure 9) also 
indicates the close correlation between the field-induced 
changes in p and m. Indeed, the comparison of Figures 4 
and 9 shows that the magnitude of resistivity correlates 
with the change in the magnetization: the suppression 
of the resistivity is large around Tc where the largest 
growth of m is seen (see Figure 4) but is weak far below 
and above Tc- Our calculation of the H dependence of 
p is in good qualitative agreement with experiment (see 
Figure 7 in Ref. |l|). 

We define the magnitude of the relative magnetoresis- 
tance as 

Sp ^ p{T,0)-p{T,H) 

p p(r,o) ■ ^ 

In this definition, 6p/p is positive and cannot exceed 
100%. Figure 10 shows the magnetoresistance as (a) a 
function of H and (b) as a function of for differ- 
ent disorder strengths at T = Tj,- It is seen from the 
comparison of Figure 3 and Figure 10(a) that the field- 
dependence of the magnetization is directly reflected in 
the field-dependence of the magnetoresistance. The A 
dependence of 6p/p at fixed H is shown in the inset to 
the Figure 10(a). The comparison of this inset with the 
inset to Figure 3 shows that disorder dependences oiSp/p 
and m at fixed H are also similar. In the weak-disorder 
regime, Sp/p is slightly decreased. Although the double- 
exchange mechanism of ferromagnetic ordering becomes 
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weaker in this regime, it is still active. In the moderate- 
disorder regime, strong-coupling-induced ferromagnetism 
(discussed in Section IV) replaces the double exchange, 
and the magnetization begins to increase with increasing 
disorder. This leads to an increase of the magnetoresis- 
tance for both the moderate-disorder and strong-disorder 
regimes. At A = 1.0, 5p/ p can attain values near 70% in 
a weak field of = 0.01. 
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FIG. 10. (a) Magnetoresistance 5p/ p as a function of the 
external magnetic field H and (b) magnetoresistance as a 
function of for difi^erent disorder strengths. The inset 
shows the disorder-dependence of 5p/ p for H = 0.005. The 
line styles in (a) are the same as in (b). 

Figure 10(b) shows that the magnetoresistance can be 
expressed by a scaling law 



— — Cm , 



(46) 



where the scaling constant C is independent of T only 
for the pure double-exchange system (A = 0). The re- 
lation ( ^ ) is approximately satisfied for finite disorder 
< A < 1.0 (at least for < 0.01), but the coeffi- 
cient C is rather high, about 4, for A — 1.0 (n — 0.8). 
versus its value of 1.9 for weak-to-moderate-disprder 
(0 < A < 0.6). Note that Furukawa's calculationcZl per- 
formed for the pure-double exchange system with clas- 
sical local spins gives the value of 4 for C (when n is 
also equal to 0.8). The difference between our estimate 
of C and Furukawa's arises from the different density of 
states. Note that Kubo and OhataQ obtain C = 1 for the 
quantum double-exchange system. The scaling constant 
C also depends on band filling. C decreases with decreas- 
ing electron filling. In particular, C = 1.7 (A = 0) and 
C = 2.5 (A = 1.0) for n = 0.67. Thus, we can conclude 
that the band filling near n = 0.8 and relatively strong 
disorder (A > 0.9) are the most favorable conditions for 
a large magnetoresistance. 



VII. CONCLUSION 

We have considered the influence of diagonal disorder 
on a simple double-exchange model with local Ising spins 
by employing dynamical mean-field theory. We choose 
a binary probability distribution for the disorder which 
greatly simplified the analysis and allowed us to exam- 
ine directly the MIT. The manganites are too compli- 
cated a system to be described completely by this sim- 
ple model. Nevertheless, we still arrive at some useful 
conclusions: (i) double exchange alone cannot explain 
the metal-insulator transition in manganites (in the best 
case, it can be applicable to LSMO at x ~ 0.3 which has 
a relatively high Curie temperature and displays bad- 
metal behavior in the paramagnetic phase); (ii) double 
exchange plus disorder cannot explain the temperature 
dependence of the thermopower S{T) because it does 
not yield a change of sign in the paramagnetic phase. 
An explanation of the large peak in S(T) (in the ferro- 
magnetic phase) is beyond dynamical meau. field theory, 
because it does not include magnon dragcS; (iii) the ef- 
fect of diagonal disorder (induced by chemical substitu- 
tion) is important, and can infiuence the properties of 
the material if the disorder strength is large enough to 
be close to the MIT; and (iv) we identified three disorder 
regimes which display different characteristic behavior. 
The weak-disorder regime has little effect on the system 
and just renormalizes the double-exchange mechanism of 
ferromagnetic ordering. The moderate-disorder regime, 
corresponds to the transition region between weak and 
strong disorder. The interacting density of states devel- 
ops a pseudogap that promotes behavior similar to ther- 
mal activation. The double-exchange mechanism is grad- 
ually replaced by a strong-coupling ferromagnetism (see 
Section IV) which has a mean-field-like magnetism. In 
this regime, the Curie temperature is sharply decreased 
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and the temperature dependence of the resistivity reveals 
a MIT at r = Tc- The strong-disorder regime, is char- 
acterized by a gap in the interacting density of states 
at high temperature and the residual resistivity exceeds 
the resistivity at T = Tc- Ferromagnetic ordering of lo- 
cal spins (with small values of Tc ^ 1/^) occurs due to 
strong-coupling physics. This regime is most favorable 
to obtain a large magnetoresistance. 
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